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Abstract 



Real dusty crystals are inhomogeneous due to the presence of external forces. 



a 

(\ [ We suggest approximations for calculations of different types of inhomoge- 

j>^, neous DC (chain and DC with a few slabs) in the equilibrium state. The 

^ OhI results are in a good agreement with experimental results and can be used as 

an effective diagnostic method for many dusty systems. 
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Formation of dust crystals (DC) takes place in a vertical electric field of the sheath, the 
gravitational field and a horizontal electrical field. The external field, acting in vertical and 
horizontal traps, stabilizes the S-dimensional DC of finite size and linear chains of d-ions 



^ ' (horizontal traps for confinement of one- dimensional DC are used in [0,0). The pressure 



of the boundaries and the external field violate the translational invariance and lead to a 
dependence of the distances between nearest neighbors in the lattice of dust particles on the 
position of the particles (see Fig.|I] and Fig.0). Therefore macroscopic inhomogeneity in a 
lattice is a new phenomenon not present in the usual infinite (very large) crystal. 

Even in the approximation of central force for interparticle interaction between d-ions, 
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DC possess a layered structure (the layered structure of usual atomic crystals, as graphite, 
is connected with the anisotropy of the interparticle interaction). 

The vertical and horizontal distances between nearest neighbours (lattice "constants" R\\ 
and R±) are in general different functions of position in different directions from the center 
of the crystal (center of inertia). Deformation of DC in the fields of the traps depends on 
its characteristics and on the plasma parameters. Therefore the electrostrictional response 
of d-ion systems on a static external disturbance can be used as a diagnostic tool for DC 
and the surrounding plasma. In particular the charge Q of tZ-ions, the screening length Ro, 
the concentration of the small ions and the electric field in the sheath can be determined. 
In the present paper the possibility to use the inhomogeneity of DC for plasma diagnostics 
is considered theoretically. 

Recently dusty plasma diagnostics appear on basis of investigations of the dispersion 
curves uj{k) for d-ion sound and properties of forced oscillations of linear (i-ion's chains in 
an electric field [|I| and under the action of laser impulses . The static diagnostic, suggested 
in this paper is simpler for the theoretical description and experimental realization than the 
dynamic sounding considered in |]I]-§]. 

For the description of a lattice configuration of A^ (i-ions in a state of deformation under 
action of external gravitational and electric forces /„ = —Wn and interparticles forces 
Fn = —VUn we will use the balance equations. Here Vn is the potential energy of the d-ion 
with number n, Un is the potential energy of interaction between the rf-ion with number 
n and all other ones. We do not take into account the force connected with momentum 
transfer from the small ions to the tZ-ions. This force very often can be omitted, because in 
the case when it is essential, the d-ions can be found not only below the sheath, but also on 
top of it, which is not observed in the experiment discussed below. 

We will use the simple approximation of nearest neighbours for the description of inter- 
particle interaction. This approximation apparently gives a good picture of the inhomogene- 
ity of DC under the action of external forces and with a screening length Rd ~ i?^, R\\. We 
also will neglect a possible dependence of the d-ion charge Q on the location in the inhomo- 



geneous DC (w.r.t. d-ion density). Therefore we suggest Q = const in our considerations. 

II. EQUATIONS OF STATIC EQUILIBRIUM 

For the case of inhomogeneous three-dimensional DC we will use a simple quasi-one- 
dimensional model of DC, in which the layer lattice with a real potential is changed into 
a one-dimensional vertical chain of particles. The effective potential for this model can be 
calculated by integration of the interaction with the nearest layer with distributed charge 
a = Q/Sq {Sq is the surface for a d-ion in horizontal direction) 

oo 

(f/(r))., = ^ I dppu (v/i^T7) (1) 



For the Debye-Hueckel interaction and simple hexagonal lattice the potential (|ip has the 
form 

This model permits to calculate the dependence of the distances between the nearest 
slabs Rn{z) as a function of height. 

In the general case of pair interaction between the d-ions in the external electric and 
gravitational fields of the sheath the potential energy can be written in the form 

Af-l N 

U + V=J2Uk + J2^''^ Uk = U{Rk), Vk = V{zk), Rk = Zk+i-Zk (3) 

k=l k=l 

Here we take into account only interaction between neighbouring particles. The potential 
energy for a horizontal chain of A^ interacting (i-ions in the electric field of the trap, has an 
analogous form and stabilizes this chain in the x-direction {zk -^ Xk). The conditions of 
balance of external and internal forces lead to a system of equations which determines the 
configuration of the rf-ions: 

UL-U',^, + VU, = 0, fc = l,2,3,...iV-2, U', = ^^ ^^' = ^' 

-U[ + VI = 0, (4) 

f/;,_i + v;^ = 0. 



Summation of the left parts of these equations leads to the obvious condition of zero sum 

N 

of the external fields: Yl ^k — 0- 

fc=i 

For the stabilization of horizontal chains an external field in the form of a parabolic well 
in the chain direction has been used in ^^. 
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14 = -mt^o(x„ - Xo)^ ^0 = j;^^Xk (5) 

fc=i 
Here Xq is the center of inertia for a chain and uq is a parameter, which determines 

the shape of the pit. According to |0] the vertical electric field in a sheath changes linearly 

with the height. This dependence is realized approximately in the regions not too close 

to the lower electrode and the border of the presheath: the quadratic approximation for 

the potential ip{z) in the plasma layer is also used in [^,|| for the analysis of the equations 

of motion of DC . Therefore in the case of a vertical potential well we use in eq.(4) the 

expansion 

V{zk) = mgzk + QifQ + Qv5o(-2fc - Xq) + \ujl{zk - X^f, 

(6) 
9^0 = V5(^o), V^o = V5'(^o), mujl = Qip"{Xo). 

The parabolic approximation for the vertical electric field is reasonable for the case of 
sufficiently thin DC . To estimate the maximal thickness i = z^ — Zi = 2{Xq — zi) of DC, 
for which this approximation is true, let us consider ip{z) = ip{0) exp^—z/Rn) and use the 
condition 

^ W^l (Xo - z,f ~ i W,\ (^) < Q^'^iXo - z,) ~ W,\ ^ 

Then the necessary inequality is £ < QRd which is usually satisfied (see, for example, 
0J^J^). The linear terms in Eq.(|^) are really absent because of the condition of zero total 
external forces: 

mg + Qip\Xo) = Q. (7) 

This condition determines the position of the center of inertia for the system of levitated 
(i-ions. 



By use of the parabolic approximation (^) in the balance equations (^) and subtracting 
from each equation the previous one, we find 

2f/fc - f/fe+i - UU + muolRk = 0, A; = 2, 3, ...N - 2, 
< 2U[ - f/2 + mulRi =0, (8) 

2f/^_i - f/^_2 + mtu^i^^v^i = 0. 

As follows from eq. (H) the intervals Rk are symmetric with respect to the center: 



Ri — Rn-i, R2 — Rn-2i ■■■■, Rk — R 



N-k, 



III. STRUCTURE OF DC WITH AN ATTRACTIVE (FOR LARGE DISTANCES) 
AND WITH A PURELY REPULSIVE POTENTIAL 

According to eqs.(^), (H) for isolated systems of (i-ions (Vl = 0), there are two different 
possibilities when external fields are absent. 

If the pair interaction between (i-ions is a nonmonotonic function and leads to repulsion 
at small distances and to attraction at large distances, then the solution of eq. (^ reads 

u[ = u^ = ... = f/;_i = (9) 

This solution describes a homogeneous chain of A^ rf-ions with equal distances between 
nearest neighbours Ri = R2 = ■■■ = Rn-i = Ro- The potential energy has a minimum 
for this state. In this case the weakly inhomogeneous configurations of (i-ions with external 
force VI 7^ can be described on basis of small deformations \Ro — Rk\ <^ -Ro- Then we use 
the expansion 

U{Ru) = Uo + -^{Rk - Ro?, mn^ = U"iRo) (10) 

If the pair interaction U{Rk) has a monotonic purely repulsive form, the (i-ions of an 
isolated system are unstable and, according to (P) all Rk are infinite. In this case stabilization 



of the system in a weak external field playing the role of a trap, leads also to a slightly 
inhomogeneous state, in which the deviations of the intervals from the average are small, 

Ro = — — -^-Rfc, \Ro — Rk\ '^ Ro (11) 

fc=i 

In this case the alternative quadratic expansion of the energy for d — d interactions has 
the form 

N-l N-1 N-1 



J2Uk = {N- l)Uo + K J2(Rk - Ro) + |mfi' J2^Rk - RoY = 

k=l k=l k=l 

[N - l)U{Ro) - (iV - l)U[,Ro + U'^{zN - z,) + imfi^ ^(^^ _ ^^)2^ (12) 

fc=i 

Uo = UiRo), U;, = ^^, rnn^ = ^^^, Y.R, = z,-z,. 

d{Ro) dR^ f^ 



k=l 

For a potential with a well U'{Rq) = the expansions (10) and ([T2|) coincide, therefore 
small deformations Sk = Ro ~ Rk of the system in an external field can then be described 

by the general equations of force balance: 

2 9 

2cosht ■ Sfc - Sfc+1 - Sfe_i - — i?o = 0, cosht = l + — — , A; = 2, 3, ..., iV - 2 

2cosht-5i-.2-^i?o-^ = 0, (13) 

2 cosh t ■ SN-1 - sn-2 - 7^-Ro ^ = 0. 

Here for purely repulsive interaction Rq is the average. For the case with attraction 
Uq = and -Rq is the equilibrium distance in the isolated system of d-ions. 



IV. SOLUTIONS AND NUMERICAL RESULTS 



A general solution of the equations in finite differences (|T3|) can be obtained in the form 

Sk = RQ-Rk = Ro + Ae''' + Be-'''. (14) 

Taking into account the symmetry of the system Sk = SN-k, the connection between 
the coefficients B = Ae^' can be found. The coefficient A can be found from the boundary 
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condition for fc = 1 (or for k = N — 1). Finally for the interval number k and purely repulsive 
potential we find 

J, -IT, f^:(^^£°!Mi-^_/p U'iRo) \ CLi(cosht) + C;,.,,i(cosht) 
R,-[R. ^ I ^^r:— -[Ra- -^^^^J C'^_^{cosht) 

(15) 



mfi2 J cosh^ V mfis y C'^_^{cosht) 



Here C'^{x) are the Gegenbauer polynomials. For the case of interaction with attraction 
U'{Rq) = 0, the intervals Rk have the form: 

cosh (f - A;) t 

Rk = Ro rivT — • (16) 

cosh yt 

Therefore in the parabolic trap formed by the external forces, a chain of d -ions is 
compressed symmetrically w.r.t. the center of inertia, and the central regions more strongly 
than the ones outwards Ri > R2 > ... For the resulting electrostrictional reduction of the 
length i of an entire chain it follows from Eq. (0) that {Rq is the equilibrium distance in a 
homogeneous chain) 

Ecosh I sinh ^^^t , , , , 

R, = 2Ro—^ -V < (^ - 1)^0 (17) 

, , smh t ■ cosh ^ 

For sufficiently long (A^ ^ 1) horizontal chains and for (in vertical direction) quasi-one- 
dimensional dusty crystals the profile distributions of charge density and mass and thereby 
the "constants" of the elastic forces can be obtained in the approximation of continuous 
media by use of eqs.(|^), (|1^). The surface density of charge is proportional to to the mass 
density and therefore there is balance of the external volume electric and gravitational forces 
in each point of a horizontal plane at fixed height. This means that even inhomogeneous 
planes (Fig.[^) and horizontal chains (Fig.|^), which are more dense in the center, are not 
suspended in the central part of the dusty system, where the density is higher. Enlargement 
of the density in the center of horizontal crystalline planes is observed in the experiments [0 , 
but quantitative measurements are unknown to us. Parallel to oscillation and wave measure- 
ments in horizontal chains, the equilibrium positions of (i-ions have also been determined in 
the electric field of a horizontal trap |ll],0,@]. According to the data of these papers for the 



case A^ = 12 the ratios of the intervals between neighbouring (i-ions in the direction of the 
center are Ri : R2 : ...-Re = 1-44 : 1.22 : 1.11 : 1.05 : 1.01 : 1.00. These resuhs are reasonably 



described by our formula (p!5D, in which for t = 0.18 (and correspondingly uq = 0.2il) these 
ratios are 1.43 : 1.27 : 1.15 : 1.07 : 1.02 : 1.00. 

The experimental data for the other half of the chain R^ : R^ : ...Rn = 1 : 1.01 : 1.01 : 
1.08 : 1.20 : 1.32 agree less with our theory for the (with respect to the center) symmetric 
chain and they are essentially different from the experimental data for the first half of the 
chain. We think that this asymmetry is a consequence of the asymmetric and not exactly 
parabolic V{x) ~ ^itiuJqx'^ shape of the external electric field (here x is the distance from 
the center of the chain). According to 0] mujQ = 2.55 ■ 10~^^ kg-s~^, m = 6.73 • 10^^^ kg. 
Using the data on the equilibrium configuration _R„ and the parameter of the trap Uq = 6.15 
s~^ we find the important characteristic oi d — d interaction Q = 5ujq = 30.7 s^^. 

For the chain with A^ = 4 the experimental data, according to |[l|,§], give ujq = 6.25 s~^ 
and i?i = 1989 /im, R2 = 1910 /im, R^ = 2031 yum. The average interval R_i = 1960 yum for 
the case A^ = 4 is twice as large as i?_L = 10^ /im for the chain with A^ = 12. This is probably 
connected with the higher charges (almost three times) of the (i-ions in |]I| and therefore with 
the stronger repulsion between them at the partially same compressing external field of the 
horizontal trap. For the conditions of the experiments lU the asymmetry of the external 
field, connected with the nonquadratic form of the potential V{x) is still stronger than 
in and this was the reason to use for the bordering intervals the expression (Ri) = 
{Ri + i?3)/2 = 2010 yum. Then according to (^) we have (Ri) / R^ = I + ujll2Vt^ = 1.05 
and n = 3.IQUJ0 = 19.7 s'K 

It is necessary to emphasize that all the results for the case A^ = 4 and A^ = 12 are 
applicable for both cases: purely repulsive d — d interaction and d — d interaction with an 
attractive part, because, as follows from eqs.(|T5[) and (p^, the ratios of intervals Rk are the 
same in these cases. 

The known experimental data on equilibrium intervals -R|| between neighbouring ions in 
vertical traps concern only dust crystals with two horizontal crystalline planes (A^ = 2) and 
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u^ + 29? 



have been obtained in [^,||]. In [Q a dust crystal with A^ = 3 has been investigated but the 
thickness of the crystal was not measured. 

According to ||,^ the ratio (i?o — -R||)/-Ro = 0.2 and does not depend on the ion's mass. 

In IP experiments are reported with dust crystals, formed by rf-ions with radii 4.7 /xm 
and 2.4 /x which leads to a difference of gravitational force proportional to mi/m2 — 8. 
The position of the center of inertia Xq of the dust crystal must be considerably changed 
in this case: a lighter crystal will shift over a distance ~ Rd, as follows from Eq.(0). A 
measurement of this effect was not reported in P|. 

According to (|16D 

D , ,2 

0.2 (18) 

and therefore f2 = \/2ujq. In contrast with [l|J^ the parameter lUq = —Vq{x) for the 
vertical electric field in a sheath is here unknown. It can be determined only on basis of 
knowledge of the interaction potential between the d-ions, via the parameter fl^ = —U"{R\\). 

Purely repulsive interaction (|l|) leads to another result. For N = 2 the exact system of 
balance eqs.(^ has the form 

-f/'(i?ll) + \/'(zi) = 0, R^\ = Z2-z^, Xo = ^^^, 

^ (19) 

f/'(i?||) + F'(^2) = 0, ;22,i=Xo±ii?||. 

In this case a more general model for the external potential (|^) than the linear one can 
be used for the description of the electric field in a sheath. Let us take 

E(Xo±^i?,l)=E(Xo)e±^'^«ii. (20) 

For the position of the center of inertia we have 

mg = QE{Xo) cosh -xi?,,, (21) 

and according to ( p!9D with U{R\\) taken from (|) we find 

1 ,47rQ2 



V'{z2) - V'{zi) = 2QE{Xo) sinh -xi?,, = 2^=^e-"^ii . (22) 



2 " V3R 
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By eliminating E{Xo) from these equations we finally find 

tanli^ = ae-^ii, a= J . (23) 

Using the experimental data for m, Q, i?x = 450 yum, _R|| = 360 yum we estimate 
the Debye radius Rd = 973 yum. In the case of a dust crystal with a lower m, Q and 
i?_L = 350 yum, R\\ = 280 yum (see also p) we find Rd = 933 yum. Therefore the Debye 
length Rd is approximately of the order of the interval between nearest (i-ions, which is in 
agreement with the estimates of |P. 

For the electric field in a sheath we find from eq. (^T]) E{Xq) = 2.84 ■ 10^ V-m^^ and 
E{Xo + A) = 1.99 ■ 10^ V-m~^ for the case of d-ions with radii 4.7 /xm and 2.4 yum. 

Let us neglect small changes of Debye radius Ro and take 

Then we obtain for the shift upwards A of the lighter crystal 

A = 0.35 (Rd) = 332 yum. (25) 

The moving of a dust crystal inside the sheath can be observed by different microgravity 
experiments (see some discussion for example in 0). 

We suggest here some experiments in which the properties of DC can be studied under 
conditions of microgravity and even changing gravity. 

One of these experiments (under terrestrial conditions) can be performed in a horizontal 
discharge, where in the horizontal direction there is only an electric force and momentum 
transfer from the small ions to the dust particles. For such experiment the latter force can 
be very essential in contrast to the conditions considered in this paper. 

The second group of experiments is connected with the effective gravity created in space 
stations by rotation of dusty plasma. If h and Qes are the distance from the axis of rotation 
to the negative electrode and the acceleration of the center of inertia for the dusty system 
respectively, the obvious connection is given by 
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ges = uj\h-Xo). (26) 

For Qcs = g and h = 1 m (rotation of the container inside the space station or rocket) 
or /i = 10 m (rotation of the space station as a whole) we find uj = 3 s^^ and uj = 1 s^^ 
respectively, which are conditions of weakly inhomogeneous {h ^ Rk) artificial gravitational 
field where our results, obtained above, are applicable. Measuring the dependence Xq = 
Xq{uj) would permit to investigate the profile of the electric field in a sheath and other 
characteristics of the dusty system and plasma. Of special interest is the investigation of 
the deformation of DC in an essentially inhomogeneous rotation field {h — Xq ~ 0.05 m and 
cu ~ 15 s^^). A detailed consideration of such experiment will be given in a separate paper. 

In the case of a dust crystal with three horizontal crystalline planes the static equilibrium 
is described by the system of eqs.(^) with A^ = 3. In the approximation for the electric fields 
used before the coordinate of the average li-ions Z2 and the value Eq can be eliminated on 
basis of balance of external fields: 

3 

3mg = QEo J^ e""'" = Q^oe"""Hl + e"^^ + e~"^^). (27) 

For the system of equations determining the vertical intervals -Ri and R2, we obtain 

1 - e-""^2 



^ e 



-xiJi _ 2e-^^2) + — = 0, 



^ + lll+_r''' ' (28) 



f (e-^^2 _ 2e-^^i) + — = 0. 



Even for the highest pressure of neutrals in 0, p = 300 mTorr {Q = 7.2 ■ lO'^e, R± = 
0.28 mm, x_Rx = 0.61) the parameter a/3 = 0.053 ^ 1. Suggesting Ri = R2 = R\\ and 
xi?|| ^ 1 it follows from eq. (|28|) that 

ry 

xR\\ ^ = 0.14. (29) 

" 1+a ^ ^ 

Let us emphasize that the vertical compression is symmetric (-Ri = R2) with respect to 
the central plane only for xR^ ^ 1. In contrast to the approximate equations (^ for the 
parabolic wells, the exact equations (53) are not symmetric for the interchange Ri ^ R2- 
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From eq.(p7|)-(p9[) with the Debye radius given in it follows that vertical compression 
is important: 

R I — ixii , . 

-^^ = 0.77. (30) 

In the framework of the quadratic approximation for the potential energy of the system 
with A^ = 3 we find according to eqs. (p!5|) -(|l6|) 

R^ ui + n^- ^ ^ 

Unfortunately the vertical interval R\\ has not been measured in |0|. The experimental 
data obtained in are not sufficient to choose which variant is preferable for purely repulsive 
interaction or interaction with an attractive part: the quadratic model or the more exact 



description (l^)-®. 



V. CONCLUSIONS. 

The method of dusty plasma diagnostics discussed above and based on an analysis of the 
inhomogeneity of the linear structures of tZ-ions, seems very attractive. In contrast to the 
situation in the usual sound method, the (i-ions of a small dust crystal or a linear dust chain 
have additional degrees of freedom. It gives the possibility to extract additional information 
from the static response (change of the equilibrium distances between the d-ions) or the 
dynamic response (oscillations and waves in inhomogeneous structures). The sounding by 
small clusters of (i-ions cannot change essentially the plasma parameters (although some 
distortion of the micro-field in the plasma can be stimulated by the traps, which stabilize the 
(i-clusters) . The advantage of static diagnostics is the simplicity of the measurements of the 
inhomogeneous structure and the simple connection with the parameters of the interaction 
between d-ions, their shielding and the characteristics of rf plasma. The precise theoretical 
consideration of the dynamical experiments [|T],^ , which are based on the excitation of the 
eigenmodes in linear chains and dust crystals, seems a more complicated problem. 
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We would like to stress that the most general consideration of the equilibrium inho- 
mogeneous configurations of dusty systems can be based on translationally non-invariant 
solutions of the connected system of kinetic equations for plasmas and Poisson's equation, 
where the separation between an external field and d — d interaction is absent. The equilib- 
rium positions for the tZ-ions can be found as the points of space where the self-consistent 
electric field is in balance with gravity. However, this program is too complicated and, as we 
showed, not necessary for a reasonable theoretical description of the existing experiments. 

We are grateful to Dr. A.Melzer for fruitful discussions and private communications 
about the experimental results. We also would like to thank Dr. H.Thomas and Dr. J.Goree 
for private communications connected with the papers [pHTT[|. 

This work have been performed with the support of INTAS grant N 96-0617. 
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FIGURES 

FIG. 1. An example of inhomogeneous plane dusty crystals (2D crystal). 

FIG. 2. Linear chain of dust particles (number of particles A'^ = 12) in a parabolic trap. The 
inhomogeneity was calculated on basis of Eq. (15) for t = 0.18 and compared with the experimental 
data [1]. 
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